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PRELIMINARY SUBMILLIMETER SPECTROSCOPIC MEASUREMENTS 
USING A SUBMILLIMETER HETERODYNE RADIOMETER 

ABSTRACT 

A submillimeter heterodyne radiometer has been under development at the Goddard Space Flight 
Center for several years. The radiometer uses a submillimeter laser, pumped by a COj laser, as a 
local oscillator and a room temperature Schottky barrier diode as the first IF mixer. The radiometer 
can resolve spectral lines in the submillimeter region of the spectrum (arising from pure rotational 
molecular transitions) to within 0.3 MHz, using a newly developed acousto-optic spectrum analyzer 
which measures the power spectrum by simultaneously sampling 0.3 MHz wide channels over a 
100 MHz bandwidth spanning the line. This report describes preliminary observations of eight 
spectral lines of H 2 O 2 , CO, NH 3 and H 2 O, all lying in the 434-524 micrometer wavelength range; 
all eight lines were observed using two local oscillator frequencies obtained by operating the sub- 
millimeter laser with either methyl fluoride (CH 3 F) or formic acid (HCOOH) as the lasing gas. 
Sample calculations of line parameters from the observed •.!. . a show good agreement with estab- 
lished values. One of our development goals is the size and weight reduction of the package to 
make it suitable for balloon or shuttle experiments to detect trace gases in the upper atmosphere. 
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INTRODUCTION 


For the past several years an effort has been underway in tlie Instrument Electro-optics Branch 
of the Goddard Space Flight Center to develop a submillimeter heterodyne radiometer using a sub- 
millimeter laser as the local oscillator and a room temperature Schottky barrier diode as the first 
IF mixer. 

The origins and development of this project are described in Reference 1 ; in this report we will 
give just a brief description of the experimental apparatus and an overview of the spectroscopic 
results obtained, with some analysis of the more significant spectroscopic measurements. 

During the past year the radiometer has been used to make measurements in the laboratory 
of selected pure rotational lines of several molecules. These spectroscopic investigations have a 
three-fold purpose: 

1. The main purpose is to demonstrate that the radiometer functions as intended - that it is 
capable of producing high resolution spectral information with reasonable integration times. 
The accuracy of the spectral data is assessed by comparing it with independently known 
results. 

2. The second purpose involves astronomical investigations of planetary atmospheres and 
the composition of interstellar nebulas. Such work is currently underway at Goddard; this 
instrument has been used with the NASA three meter infrared telescope at Mauna Kea, 

Hawaii for studies of tlie Orion nebula, under the direction of scientists from the Infrared 
and Radio Astronomy Branch. The spectral data obtained in the laboratory are useful in 
interpreting the results of these field experiments. 

3. A third objective of the investigations is the support of upper atmosphere research, 
especially studies of the formation and breakdown of ozone. The radiometer should be able 
to detect trace amounts of various gases involved in these processes. Laboratory spectro- 
scopic data are needed for both the design and analysis of such experiments, because 
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the atmospheric transmission problem must be solved for both tlie simulation of observations 
(to determine detectability) and the inversion of actual observations (to retrieve actual con- 
stituent concentrations). Spectroscopic data at submillimeter wavelengtiis do not generally 
exist; instruments of the type described in tliis report are tlie only e.xisting means of acquiring 
tlie needed lugli resolution data. 

Preliminary measurements have been made on selected spectral lines of HjOa, CO, NH3 and 
HjO; these are listed in Table I, 


Table 1 

Spectroscopic Lines Investigated 


Molecule 

Line-Center 

Frequency 

(MHz) 

Quantum 

Transition 

! 

Pressure 

Range 

(mtorr) 

Receiver 

Noise 

Temperature 

/V\ 

VW 

HaOa 

599,723.63 ±0.1 

9(0,9) ^9(1,9) 

300 

i 


HaOa 

601,885.28 ±0.1 

11(0,1 D'f- 10(1,9) 

5-1,800 

5,000 

HaOa 

603,394.93 ± O.I 

8(1,8) ^8(0,8) 

50-2,400 

6,000 

HaOa 

606,717.54±0.1 

7(1 ,7) ^7(0,7) 

5- 1,800 

5,000 

HjOa 

608,865. 

21(3,19) 22(2,21) 

300- 1,800 

8,000 

CO 

691,472.97 ±0.1 

J = 5 -^6 

2- 1,600 

10,000 

NH3 

572,498.15 ±0.15 


5-650 

20,000 

HaO 

620,700.81 ±0.39 


1,000 

25,000 


These lines all arise from pure rotational or rotational-torsional transitions in the vibrational 
ground state. The line-center frequencies are taken from Reference 2, with tlie exception of the 
608865 MHz line of H3O2 , for which the line-center frequency given above was measured in our 
laboratory. The frequencies measured in our laboratory for the otlier lines agree with the values in 
Table 1 to within a few MHz in most cases; tliis is quite good agreement, considering that frequency 
measurement was not a primary aim in our experiments. For the 608865 MHz line, the value given 











in Reference 2 was calculated from theory ; it differed from our measured value by several hundred 
MHz. 

For each of the lines studied, the nearest neighboring lines (according to Reference 2) are 
separated from the line investigated by mors than 1,000 MHz. Since the largest bandwidth of any 
of the spectrum analyzers used is about 300 MHz, it is clear that each of the lines studied is effec- 
tively isolated; i.e., not overlapped by any neighboring lines. 

The range of sample cell pressures used allowed observations of the lineshape in both the 
doppler and pressure-broadened regimes, for some of the lines. 

From the preliminary data obtained with the heterodyne radiometer and from comparisons of 
some of the data with theoretical results it was found that: 

1. The radiometer is capable of actually achieving resolution to better than the doppler 
width; the best resolution to date is 300 KHz. 

2. The radiometer is capable of giving the kind of laboratory spectroscopic data required 
for the design and analysis of upper atmosphere experiments. 

3. Comparison of the radiometer data with some theoretical results and with previously 
known data showed good agreement. 
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EXPERIMENTAL APPARATUS 

The absorption line to be observed is produced by directing the radiation from a hot black- 
body source through a long cylindrical sample cell filled with the absorbing gas at low pressure. It 
is also possible to produce an emission line by replacing the hot source by a cold background; for 
example, an absorbent material cooled in liquid nitrogen. 

The first heterodyne stage is a quasi-optical arrangement which translates the spectral line(s) 
bemg observed from the submillimeter region of the spectrum down to the microwave region. 
The local oscillator is a submillimeter laser, optically pumped by a COj laser. The nearly mono- 
chromatic output of the submillimeter laser is optically combined in a diplexer with the beam 
emerging from the sample cell. The combined beam is then focused onto the whisker antenna of 
The Schottky diode mixer, which produces a down-shifted spectrum in the microwave region (18 
GHz maximum). 

The second heterodyne stage is a conventional microwave mixer which shifts the absorption 
line spectrum further down into the frequency range from 0.1 to I GHz (depending on the {yy", of 
spectrum analyzer used), where it can be precisely measured.. 

The measurement of the spectral line, now in the sub-GHz region of the spectrum, is accom- 
plished by observing it simultaneousiy in many channels, which span the width of the line. Two 
types of spectrum analyzers were used for this purpose: a conventional filter bank type and two 
different acoustooptic devices. 

Details on the development and characteristics of the suhmillimeter laser local oscillator, the 
COa pump laser, the diplexer, the Schottky diode mixer and related devices that make up the 
heterodyne receiver and data processing electronics are given in References 3-22. The laboratory 
apparatus is shown schematically in Figure I. 
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Figure 1. Submillimeter Wave Heterodyne Radiometer 

The optical pump source for the submillimeter laser is an Apollo Model 550A 9-1 1 micrometer 
CO 2 flow-through laser using a C02-6%:N-18%:He-76% gas mixture (see Reference 1). For these 

j 

experiments the CO 2 laser was operated at a pressure of 22 torr. The laser is presently capable of 
CW output up to 60w, but was operated at about 30w for these experiments. A cavity stabilizer 
PZT supplied by Apollo is used to maintain the pump line and under most experimental conditions 
a Burleigh etalon is used to lock on the pump wavelength to lend additional stability over long 
time periods. This is accomplished by monitoring a reflected beam from a ZnSe Brewster window 
located in the output beam of the laser. The Apollo has a diffraction grating with both horizontal 
and vertical adjustments to aid in locating and maximizing the required pump line. The CO 2 laser 
was calibrated with a spectrum analyzer to correlate the grating dial micrometer positions with 
pump lines in the 9 and 10 micrometer CO 2 branches. 

The subrndlimeter laser cavity is the culmination of several years of design and construction 
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by N. McAvoy, G. Koepf, C. Rossey and C. Peruso (see References 1, 6, 10, 1 1, 13, 14, 21). Tlie 
device is optically pumped by the Apollo COj laser and is designed to function as a local osciliator 
in the 400-800 GHz (~750-375 /im) range. A 1.76 meter submillimeter cavity is used to produce 
an effective 14.08 m path length by reflecting the COj beam on four round trips inside the cavity 
through the absorbing laser gas. The final reflection exits off-axis through the input port and is 
monitored by a pyroelectric power meter. The submillimeter laser beam is coupled through a 
12mm quartz window at the exit port. There are adjustments for cavity wavelength and for internal 
mirror alignment external to the cavity, located at the exit end of the cavity. Maximizing the power 
of the unabsorbed COj beam on the power meter is the technique used in making final alignment of 
the reflecting mirrors inside the submillimeter cavity. Bleed-in valves permit entry of a vaporized 
liquid from an external storage tube or direct entry of a gas simultaneously into the ends of the 
submillimeter cavity. Pressure in the cell is monitored near the middle of the cavity by a MKS 
baratron pressure gauge connected to a digital voltmeter (DVM). The cavity pressure ranges from 
about 30 to 150 mtorr, depending on the gas used. The design and operation of the submillimeter 
laser local oscillator and the COj pump laser is given in greater detail in Reference 1. 

A two meter glass ceil made in-house serves as the sample cell. It is designed to accept gas 
samples directly or vaporize liquid samples from an external storage tube. A blackbody radiator 
serves as the source of energy for absorption. Although the blackbody produces a 1,300 K source, 
present indications are that optical coupling causes losses tliat result in an effective temperature 
difference of 250 K with respect to room temperature, Pressure in the cell is monitored near the 
exit port by a MKS baratron pressure gauge connected to a DVM. The blackbody radiation is 
chopped prior to entry into the sample cell at a variable rate dependent on the filter bank require- 
ments. 

The beam emerging from the sample cell and the submillimeter local oscillator beam are chan- 
neled into a diplexer (References 1 , 9) which serves to combine the beams for input to the Schottky 
diode mixer. The diplexer is a Mach Zehnder interferometer, as shown in Fig. 2. A translator 
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Figure 2. Mach-Zchnder Type Beam Diplexer 

assembly on the diplexer permits adjustment of the reflective surface to tune the fringe pattern for 
constructive interference on the appropriate IF wavelength. This stage may also be used to accu- 
rately determine the wavelength of the submillimeter laser beam, since the interferometer spacing 
must also be an integral multiple of the submillimeter wavelength. 

The mixer is a Schottky barrier diode designed and fabricated at Lincoln Labs for use as a 
receiver in our frequency range (see References 1, 3, 5, 7, 8, 12, 15, 16, 17, 19, 20). It is mounted 
in a comer cube and placed at the focus of an ellipsoidal mirror. The comer cube helps to match 
the diode antenna pattern to the submillimeter beams; it is mounted on a base with three degrees of 
translational freedom. A battery-powered bias box is used to provide a 0.07 ma current to the 
diode. A connection to an oscilloscope provides a monitor point for the submillimeter signal. 
This first IF signal is then fed into the second IF. 

A low noise f)reamplifier and filter combination is used to boost the signal from the first IF. 
System noise temperatures are determined by matching of detectors to the low-noise preamplifiers. 
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Typical noise temperatures ranged from ifbout 4,000 K to about 9,000 K. For the preliminary 
observations of NH 3 and H 2 O, with extremely high first IF frequencies of about 12-16 GHz, the 
noise temperatures were in the 20,000-25,000 K range. If the first IF signal could be handled by 
the filter bank, as is the case for the conventional RF bank described later, then the signal was 
split: one channel going to the filter bank and the second channel to a crystal detector which was 
used to monitor the signal response from the sample. When one of the two acoustooptic spectrum 
analyzers was used, the first IF signal was downconverted by a second mixing stage (which uses a 
frequency synthesizer as the local oscillator) to the center frequency of the spectrum analyzer 
(either 150 MHz or 400 MHz). Double downconversion was necessary for H^O and NH 3 because 
of the extremely high first IF’s used. A D\^ and a lock-in amplifier (referenced to the sample cell 
chopper) were used to monitor the second IF output signal. This set-up provides a means of deter- 
mining the system noise temperature while running. The method of measuring the noise tempera- 
ture is an adaptation of the so-called Y-method of system temperature measurement; it is described 
in detail in Reference 1 . 

Three different spectrum analyzers were used dirnng these experiments: one conventional type 
using an RF filter bank and two acoustooptic an?’’ ■ zers. The conventional spectrum analyzer is 
described in Reference 23. It incorporates a second IF containing two tunable local oscillators, 
which together can accept RF signals up to several GHz and mix them up or down to the 1.170 
GHz center frequency of the filter bank. This analyzer is designed to give 5 MHz resolution, using ' 
64 channels. A block diagram is given in Figure 3. 

The two acoustooptic spectrum analyzers (AOS’s) both work on the same principle; functional 
diagrams (including the 2nd IF stage) are given in Figure 4. In both AOS’s the spectral analysis 
is accomplished by using a Bragg cell. In this device high frequency sound waves, created by | 

t 

coupling the 2nd IF signal into a crystal by means of a piezoelectric transducer, create within 
the crystal a periodically varying index of refraction along the direction of propagation of the 
sound waves, caused by the density variations. This grating-like pattern causes an incident 
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Figure 4. Acoustooptic Spectrum Analyzers (AOS I, IF) 








monochromatic light beam (from a HeNe laser) to bo deflected tlirough a small angle, the amount 
of deflection depending on the frequency of tlie 2nd IF signal. Since the signal contains a (band- 
limited) continuum of frequencies, the laser beam is spread out into a fan-shaped beam whose 
intensity varies with the deflection angle; in fact, the intensity is proportional to the amplitude of 
the corresponding signal frequency, so that the beam actually forms the spectrum of the 2nd IF 
signal. A lens then focuses the beam onto a linear array of integrating photodiode detectors; the 
output from these yields the signal power spectrum, with a finite resolution depending on the 
number of photodiodes. (The attainable resolution, and thus the number of diodes used, depends 
on the size of the diffraction-limited spots focused onto the diodes and on the angular spread of 
the deflected beam.) The tlieory and design of acoustooptic devices are given in References 24-28. 

The two AOS’s differ in the material used for the crystal (flint glass and lithium niobate for 
AOS I and II, respectively) and in design details. Functional diagrams for the two analyzers are 
shown in Figure 4; Figures 5 and 6 sltow a functional sketch and a photograph, respectively, of 
AOS I. AOS I was developed by the Instrument Electro-Optics Branch and AOS II by the 
Infrared and Radio Astronomy Branch; a more detailed description of AOS II is given in 
Reference 29. Each of the three spectrum analyzers has its own associated data collection and 
reduction system, as indicated in the figures. 
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l-iguro 5. Sketch of AOS I 



l•'lgllre 6. Photograph of AOS I 



LABORATORY TECHNIQUE 

Hydrogen peroxide (H 3 O 3 ) and water (H^O) samples were in liquid form. The liquid was 
placed in the external storage tube on the sample cell. The sample cell was evacuated to a reading 
on the pressure gauge which did not change and was calibrated to be 0 mtoir pressure. The entry 
port on the sample cell was opened and a sample pressure greater than required was noted. The 
entry port was closed off and the evacuation of the sample down to tlie desired pressure was ob- 
served on the pressure gauge. At this point the exit port of the sample cell was closed off. For 
carbon monoxide (CO) and ammonia (NH ) cylinders of gas were hooked up to the exit port of the 
sample cell. An excess pressure was added to the sample cell, then the excess was removed while 
monitoring tlie pressure gauge. It should be noted that, as the desired pressure was approached, 
the bleed-off speed was slowed considerably to reduce any major gradient effects inside the sample 
cell after reaching the desired pressure. The run was made, then the sample cell was evacuated to 0 
mtorr and an evacuated cell run was made. This technique was used to help make the measure- 
ments independent of tlie receiver noise temperature, as described in the Appendix. 

Two submillimeter laser gases were used In these experiments: methyl fluoride from a gas 
cylinder and formic acid from a liquid storage tube. It was noted earlier that the gas line into the 
submillimeter laser is split, providing gas entry simultaneously at both ends of the cavity. When 
using methyl fluoride the evacuated cavity was brouglit to a pressure of about 40 mtorr of methyl 
fluoride, then the pressure was increased by about 50 mtorr by the addition of a buffer gas, sulfur 
hexafluoride. There are twin entry valves; each valve was connected to a gas cylinder. Liquid 
formic acid was placed in a glass storage tube which was connected to one of the entry valves. A 
pressure of about 80 mtorr was used for this laser gas. 

Depending upon the strength of the line observed, integration times of 1.5 to 3 minutes were 
used during the data collection runs. It has been noted that a signal run was made followed by an 
evacuated cell run. The software associated with the data reduction, on either data system, per« 
mitted results to be plotted with nonnalization by the evacuated cell run. .Any adjustments to the 
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system optics or electronics were made prior to a new run, not between the signal and evacuated 

cell runs. During these runs the CO 2 laser and tlie submillimeter laser remained quite stable, due in 

large measure to the stable room temperature of about 22 C. It was generally not necessary to use 

» 

the Fabry-Perot etalon lock-in servo loop to maintain COj pump amplitude and frequency stability. 
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SOME THEORETICAL BACKGROUND MATERUL NEEDED TO REDUCE THE OBSERVED 
DATA TO LINE STRENGTHS AND WIDTHS 


Relation Bet^veen the Observed Data and the Absorption Coefficient 

The first step in reducing the data is to calculate the absorption coefficient a from the meas- 
ured quantity, which we will call Q. (Q is a function of the signal and reference voltages at the out- 
put of the spectrum analyzer, as defined below.) A careful analysis of the processing of the input 
signals througli the radiometer leads to the following relation (see the Appendix for a derivation); 




(I) 


where the measured quantity Q is defined as: 

o^JSlWR 

(S«-R)/R 

and a is the absorption coefficient in cm“‘ , £ is the length of the gas sample cell in cm, e is the 
frequency deviation from the line-center frequency, R is the voltage v measured at tlie output of the 
spectrum analyzer (for a given frequency interval, the widtli of which depends on the resolution 
of the spectrum analyzer) when the detector sees a reference source at some temperature (in our 
case, usually a chopper blade at room temperature), S is the same voltage when a signal is present 
(i.e., an absorption or emission line seen by the detector through the gas-filled cell) and Sq is the 
same quantity, with the cell evacuated. To retrieve the signal (the spectral line) from the accom- 
panying noise, the quantities S, Sq and R are digitally integrated in the minicomputer associated 
with the spectrum analyzer. 


It is important to note (see Appendix) that Eq. 2 leads to Eq. 1 only if the elapsed time 
between the measurements of S and Sq is small enough so that there are no significant drifts in 
the gain or loss factors associated with the detector, or in the noise temperature of the detector, 
Note also that Q is unchanged by any linear transformation of v (i.e., v' = av + b); thus Q is un- 
changed by changes in the scale factor or zero point of the output voltages from the spectrum 
analyzer. 


After the absorption coefficient a(e) has been obtained from the measured data Q(e) with the 
help of Eq. (1), we try to obtain from a(e) what information we can on the line strength and width. 
To do this we need explicit formulas for a(e). These are collected below. 

General Formula for the Absorption Coefficient (for a Single Gas) 

The formula for the absorption coefficient for a given line has the following general form: 

a(63>,T) = S(T)*p*g(e33,T) (3) 

where S(T) is tlie line strength per unit pressure (i.e., the integrated absorption coefficient over the 
line width, per unit pressure), T is the absolute temperature, p is the pressure (assuming only a 
single gas) and g is the normalized line shape (described below). (Detailed derivations of these 
formulas may be found in References 30-32.) 

Formulas for the Lineshape 

The lineshape g(e,p,T) is normalized so that the integral of g over the frequency width of the 
line is unity. The expressions for g for the different regimes may be put in the following forms: 

For pure doppler broadening (Gaussian shape) ; 

g/«Ii.T)= (4) 

where is the half-width of the doppler-broadened line. It is given by the formula: 

Aj (T) = i (2 kNo In 2)« n, CT (5) 

where k is Boltzmann’s constant (numerical values and units are given below), Nq is 
Avogadro’s number, Vq is the line-center frequency and M is the moiecular weiglit of the 
absorbing gas. 


For pure pressure broadening (Loren tz form): 


g,(e,P,T) 


^1 

;r(e^ +A^) 


( 6 ) 
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where is the half-width of the pressure-broadened line. It may ’ expressed in the 
form: 

A,-p*Ao(T) (7) 

where Aq is the linewidth for unit pressure. Usually the calculation of Aj (or Aq ) from 
theory is very difficult and can only be done approximately; consequently it is necessary 
to resort to experimentally measured values of A,. However, it is known that Ao(T) 
varies with T as 1/T, for a long-range intermolecular force law (see Ref, 30, page 369). 


For combined doppler and pressure broadening: 






(0 


dt 


-d i tr+p^Ag) 
where t is a variable of integration with the dimensions of frequency (sec"* ). 


( 8 ) 


General Formula for the Line-Center Absorption Coefficient 

The formula for the lineshape with combined doppler and pressure broadening is not 
analytically tractable, it may however be cast in various forms suitable for numerical computa- 
tions (see, for example, Reference 31). For the special case of absorption at the line center (e = 0), 
the formula may be put in a particularly simple form in which the definite integral reduces to the 
probability integral. To derive this form we use the following definite integral evaluation: 


where Re /3 > 0 and 1 arg m I < ff/4 (see Reference 33, page 338, formula # 3.466) and is the 
probability integral: 



X 



o 


If we apply this representation to Eq. 8, with n = (ln2)^ / A^ and 0 = pAg , we obtain the follow- 
ing expression for the line-center absorption coefficient: 
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«(0;p,T) = S(T) X e*’ ( 1 - $rn» (9) 

where x, which is essentially the ratio of the collision width to the doppler width, is proportional 
to the pressure: 

X = (In 2)^ ^ = (In 2)'' (10) 

Ad Ad 

This formula holds for all x (i.e., all pressures). 


The probability integral may be calculated from the series representation 

4>(x) = 


- 1 . ( 

V 3.1! 


_ 

^ 5.2! ~ 7.3! 


( 11 ) 


which is convergent for all x but conveniently calculable only for x « 1 , or from the series 


<J>(x) 


= 1— IT 0“ 

X \ 


1 1.3 

+ 


1.3.5 


•••) 


2x" (2x^)2 (2x^)3 

which is semiconvergent for x » i. For values of x near unity, tables of may be used. 


( 12 ) 


Equation 9 may be used to calculate tlte line strength S from experimentally measured values 
of the line-center absorption coefficient a (0; Po, Tq) at pressure Po and temperature Tq. Note 
however that the pressure broadening parameter Aq must also be known — probably from experi- 
mental observations, since it is difficult to calculate accurately from theory. 

Formulas for the Line Strength 

The calculation of the line strength *’'om theory can be very difficult, depending on the mole- 
cule, and only order-of-magnitude estimates may be obtainable. The theoretical expression for the 
line strength may be put in the following general form: 

Stt^ 

SCT) = - (1 - e-'-'-'''’') f, (T) I n„ I" (13) 

3c h k T ^ 

where h is Planck’s constant, Vq is the line-center frequency, fj is the fraction of molecules in the 
lower state of the transition and Mu are matrix elements related to tlie dipole moment of the 
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transition. The difficulty lies in calculating fj and | Some general results are available, how- 
ever (see, e.g., Ref. 30). For example, we may state the following results: 

1 . For diatomic (or linear) molecules, for pure rotational transitions in the uth vibra- 
tional state and for low J-values (J is the rotational quantum number): 


‘ 2(J+1) kT “ 


(14) 


and 




(15) 


where f„ is the fraction of molecules in the uth vibrational state (u = 0 for the ground 
state) and [x is the dipole moment of the transition (see Ref. 30, Chapter 1). It 
follows that: 


f, (16) 

for these types of transitions in diatomic or linear molecules. In particular, f„ s 1 
for the ground vibrational state at ordinary temperatures, so that for transitions 
wunin the ground state 

f. (17) 

2. For symmetric-top molecules, for pure rotational transitions within the ground 
vibrational state and for low rotational quantum numbers (J and K), 

fj varies with T as l/T^/^, approximately. 

3. For transitions involving excited vibrational states, or for transitions between high- 
energy rotational levels of asymmetric- top molecules, the temperature dependence 
off I is more complex. 
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If a good theoretical estimate of the line strength is not available, we must resort to an experi- 
mental measurement. If the line strength is strongly temperature dependent this would require 
measurement over the temperature range of interest. However we may find the line strength S(T) 
as a function of temperature from only a single measurement, provided we know the form of the 
function fi (T), expressed approximately as a power of T. For example, it is easy to derive formulas 
of the, following type from Equations 3 and 13: 

Po g(o :Po .To ) T f t (To ) 1 - e"'’ • 

This equation gives the temperature dependence of the line strength, provided we have an approxi- 
mate expression for fj asapowerofT,sayT®, because for this form the unknown coefficient of T“ will 
cancel. We need only measure the line center absorption coefficient a (0;po ,Tq) at any convenient 
pressure po and temperature Tq. In addition, we must use the lineshape function g which is appro* 
priate; for example, if Po lies in the pressure-broadened regime, we must use gj as given by Eq. ( 6 ), 
evaluated at e = 0, p=Po and T-Tq. Note that Eq. IS, once evaluated, applies to all pressure 
regimes, not just the one to which g applies. For example, if po lies in the pressure-broadened 
regime, and we use gj from Eq. ( 6 ), the resulting form of Eq. 18 applies to the doppler regime also. 
Note also that we may need, in addition to a measurement of a, a measurement of A,, since the 
pressure broadening may not be accurately calculable from theory. 


Condition for the Detectability of a Line 

The minimum resolvable temperature difference of a radiometer of our type is given by the 
formula (see Reference 34, page 102): 




(Br)« 


> 


where (AT)jjif, is the minimum detectable temperature difference (in K), T 3 YS is the system 
noise temperature (the sum of the antenna noise temperature and the receiver noise temperature) 
referred to the antenna terminals (in K), B is the predetection bandwidth (in Hz) and t is the post- 
detection integration time (in sec). If we interpret AT as the apparent line strength seen by the 
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detector, i.e., as the off-line temperature minus the line-center temperature, then the above form- 
ula relates the minimum detectable apparent line strength to the integration time, if B and ^SYS 
are regarded as constant. 
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• Units and Constants 


For all of the above equations the units and constants are as listed below. 
Units 

length - cm 
mass — g 
time — s 

charge - statcoulomb 
frequency - Hz 
temperature - K 
g -- s (or ) 

a -- cm'* 

p - dyne/cm ^ 

S — s/g(orcm"* • Hz • (dyne/cm*)"*) 

fi — statcoulomb-cm 

M — g/mol 

Constants 


c = 

2.997925 X 10*“ 

cm/s 

velocity of light in vacuo 

h = 

6.626196 X 10"” 

erg-s 

Planck’s constant 

k = 

1.380622 X 10'*® 

erg/K 

Boltzmann’s constant 

No = 

6.022169 X 10” 

mole"* 

Avogadro’s Number 


Conversion Factors 
1 joule = 10’ ergs, exactly 

1 debye = 10"*® statcoulomb-cm (dipole moment) 
1 ton = 1.33322 X 10® dyne/cm’ 

1 atm = 760 ton 
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LABORATORY MEASUREMENTS 


HiOa - 603.395 GHz 

A local oscillator (LO) frequency close to 603 GHz was obtained by using methyl fluoride 
(CH3F) as the lasing gas in the submillimeter laser. It produces the 604.2973 GHz line of CH3F 
when optically pumped by the 9.55P20 line of the COj laser. This scheme of submiliimeter laser 
local oscillator operation was used for all HjOj spectroscopy reported below. Mixing of the H2O2 
line and the local oscillator frequency in the 


Schottky diode then produced a first IF signal at 
902 MHz, which was analyzed on the conven- 
tional RF filter bank. Integration times of 
about 90 seconds were used on both the sample 
run and the evacuated cell run. Data was 
collected at the following sample pressures: 50, 
100, 200, 300, 600, 900, 1200, 1500, 1800 and 
2400 millitorr. System noise temperatures were 
about 6,000 K during the runs. An example of 
the spectra obtained is shown in Figure 7, for 
a pressure of 200 millitorr. Poor calibration of 
the vertical axis precludes any quantitative 
analysis, but it is clear from comparison of this 
spectrum with those taken at other pressures 
that the line is mostly pressure-broadened at 
200 millitorr. 

H2O2 - 608.865 GHz 



Figure 7. HOOH Line at 603 GHz, for 200 
mtorr (RF filter bank) 


Mixing of the LO and sample beam produced a 
first IF frequency of 4.568 GHz, which was 
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amplified by a liquid nitrogen cooled amplifier. 
The output was sent to the conventional RF 
filter bank. Integration times of about 90 
seconds were used for both the sample and 
evacuated runs. The following sample pres* 
sures were analyzed: 300, 600, 900, 1200, 
1500 and 1800 millitorr. System noise tem- 
perature was about 8,000 K during testing. 

An example of the spectra obtained is shown 
in Figure 8, for a pressure of 300 millitorr. 
Again quantitative analysis is ruled out by 
poor vertical calibration and uncertainties in 
the pressure measurement, but it is clear from 
comparison of this spectrum with those taken 
at other pressures that the line is pressure- 
broadened at 300 mtorr. 

- 601.885 GHz, 606.717 GHz 
The first time we studied these transitions we 
based the set-up on the 601.885 GHz line. 
The first IF signal of 2.412 GHz was sent 
directly to the conventional RF filter bank. 
Because of the double sideband conversion 
and the fact that the L.O. frequency of 
604.297 GHz is nearly midway behveen the 
two HjOj lines, the separation of the two 
lines after th., first down-conversion 



Figure 8. HOOH Line at 608 GHz, for 300 
mtorr (RF filter bank) 



Figure 9. HOOH Lines at 60 1 and 606 GHz, 
for 900 mtorr (RF filter bank; lines 
not resolved) 
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(by the Schottky diode) was only about 8 MHz; these two lines were not resolved by the R.F. filter 
bank, which itself has a resolution of about 5 MHz. During this study the integration times were 
about 90 seconds. The pressures studied were: 300, 600, 900, 1200, 1500 and 1800 millitorr. The 
system noise temperature was about 5,000 K. An example of these spectra is shown in Figure 9; 
again quantitative analysis is ruled out by poor vertical calibration and uncertain pressure measure- 
ments, but again it is clear from comparison of this spectrum with those taken at other pressures 
that the lines are pressure-broadened at 900 mtorr. 

During the second study of this pair of transitions the first IF signal of 2.41 2 GHz was down- 
converted by a second mixer to the 400 MHz IF used by the acoustooptic spectrum analyzer 
(AOS II). Integration times were about 70 seconds. The following pressures were studied: 5, 10, 
20, 40, 80, 160, 320, 640 and 1 ,250 millitorr. System noise temperatures were about 5,000 K. It 
was with the AOS II resolution of 0.7 MHz that we were first able to verify the existence of both 
lines as suggested by the absorption results seen in the first study. Some of the spectra obtained are 
shown in Figures lOa-lOf. The calibration for these spectra was good; some quantitative analysis 
is given in the section on analysis of data. It is clear from visual inspection of Figures lOa-lOf 
that the two lines are mostly doppler broadened at 5 and 10 mtorr, although some pressure 
broadening is evident even at 10 millitorr. As the pressure increases the lines grow both deeper and 
broader (due to simultaneous doppler and pressure broadening) until they merge at around 100 
mtorr. From this it is clear that the lines would not have been resolved in the first study even with 
greater resolution, unless measurements were made below 100 mtorr. 
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but merged at higlier pressures) 



CO - 691.473 GHz 


ORIGINAL F;‘w'. I.*; 
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A transition of J ■ 5*-6 produces the 691.473 GHz line of carbon monoxide. A local 
oscillator frequency close to 691 GHz was obtained by using formic acid (HCOOH) as the lasing 
gas in the submillimeter laser, which produces the 692.95 14 GHz line of formic acid when optically 
pumped by the 9.27R20 line of the COj laser. This local oscillator frequency, when mixed with 
the sample beam, produces a fint IF of 1.478 GHz. Three studies were made of this CO line using 
the same basic configuration. 


The first study was done using AOS I. 
For this spectrum analyzer the fint IF had to 
be downconverted to 150 MHz. Integration 
times up to several minutes were used. Pres- 
sures of 110, 220 and 300 millitorr were in- 
vestigated. The system temperature in these 
early trials was about 9,000 K. Figure 1 1 
shows one of the spectra, taken at 1 10 milli- 
torr. These spectra were well calibrated, 
permitting some quantitative analysis; but this 
will not be given here, because the results of 
the third study are better (see the section on 
analysis of data). It is clear from comparison 
of the three spectra that the line is at least 
partially pressure-broadened at 110 mtorr. 

The second study was done using AOS II. 
For this analyzer the first IF was downcon- 
verted to 400 MHz. Integration times of 70 
seconds were used for both the sample 



Figure 1 1. CO Line at 691 GHz. for 110 mtorr 
(AOS I) 



and evacuated cell runs. Pressures studied 
were: 60, 100, 300, 1,(X)0 and 2,000 milU- 
torr. The system temperature was about 
4,000 K. Lack of calibration of the vertical 
axis precluded quantitative analysis of these 
spectra. An example is shown in Figure 1 2. 

The third study was done using AOS I, with an improved data collection and analysis system. 
Integration times of 80 seconds were used for all pressures except 2 mtorr, which required 160 
seconds. Pressures of 2, 5, 10, 20, 50, 100, 200, 400, 800 and 1600 mtorr were investigated. 

Figures 13a-13j show these spectra. The data collection and analysis system was considerably im- 
proved for these investigations (see Reference 35), and the data are quite well calibrated. In the 
section on analysis of data these results are shown to yield a value for the dipole moment of CO (for 
the transition involved) which is quite close to the value given in the literature. It is worth noting 
also that the range of pressures goes from the doppler regime (or close to it) at the low end to the 
collision-broadened regime at the high end. Note that the frequency scales in all the plots in Figure 
13 have been arbitrarily shifted to read 0 MHz near the line center. The actual frequency of the 
zero marker is noted on each plot; these frequencies were determined by prior calibration of the 
acoustooptic spectrum analyzer resolution elements with a sine wave generator. 
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Figure 1 2. CO Line at 69 1 GHz, for 300 mtorr 
(AOS II) 
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Figure 13a. CO Line at 691 GHz, for 2 mtorr 
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Figure 13b. CO Line at 691 GHz, for 5 mtorr 
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Figure 13c. CO Line at 691 GHz. for 10 mtorr 
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Figure 13e. CO Line at 691 GHz for 50 mtorr 
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Figure I3f. CO Line at 691 GHz for 100 mtorr 
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Figure 13g. CO Line at 691 GHz for 200 mtorr 
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Figure 13i. CO Line at 691 GHz for 800 mtorr 
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Figure 13j. CO Line at 691 GHz for 1600 mtorr 
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NH 3 - 572.498 GHz 

We ran initial studies on the transition 
of ammonia which produces a line at 57 2.498 
GHz. Using HCOOH pumped by the 9. 23R28 
CO] line, a local oscillator frequency of 
584.3729 GHz was produced. Mixing at the 
diode produced a tint IF of 1 1.875 GHz. Using 
a liquid nitrogen cooled preamplifier this 
signal was downconverted to the 400 MHz IF 
of AOS II. Integration times were about 180 
seconds for both the sample run and the 
evacuated cell run. Pressures of 5, 10, 20. 40, 
160 and 640 millitorr were studied. System 
noise temperature was about 20,000 K. For 
these spectra quantitative analysis was not 
an objective: the primary purpose was to • 
verify that heterodyne techniques with a 
large IF ( 12 GHz) could produce usable 
spectra. Figure 14 shows an example of the 
spectra obtained. 

H,0 - 620.701 GHz 

* 

Finally, a short study was done to see if 
the system could locate the transition of water 
which produces a 620.701 GHz line. Using 
methyl fluoride as the submillimeter laser 
gas. optically pumped by the 9.55P20 COj 



Figure 14. NH 3 Line at 572 GHz, for 5 mtorr 
(AOS II) 



line, a local oscillator frequency of 604,2973 
GHz was produced. Mixing of this frequency 
with the H 5 O line produced a first IF 
frequency of 16.403 GHz. Using a liquid 
nitrogen cooled preamp, this signal was down- 
coverted to 150 MHz by the mixers available 
at the time and then sent to the conventional 
RF filter bank, where the on-board mixer 
generated the 1.170 GHz signal it required. 
Integration times of 60-180 seconds were 
used for a pressure of 1,000 mtorr of H^O. 
Movement of the second IF allowed tenta- 
tive verification of the line, although no further 
tests were done. The system noise temperature 
was about 25,000 K. For these spectra quanti- 
tative analysis was not an objective ; the 
primary purpose was once again to verify 
that heterodyne techniques with a large IF 
(16 GHz) could produce usable spectra. A 
sample spectrum is shown in Figure 15. 



Figure 15. HjO Line at 620 GHz, for 1000 
mtorr (RF filter bank) 
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Emission Stiidies 

There were additional tests to study the emission spectra of H^Oj and CO. For these tests 
the hot blackbody source was replaced with a piece of ecosorb dipped in liquid nitrogen. These 
preliminary results support the utility of the radiometer for emission studies. An example of these 
spectra is shown in Figure 16; it is an emission spectrum of the RQO 9(0,9)'*~9( 1 ,9) line of 
at 599 GHz, taken at a pressure of 300 millitorr. 



Figure 16. Emission Line of HjOj at 599 GHz, for 300 mtorr (RF filter bank) 


34 



ANALYSIS OF DATA 


The analysis of the data given below is intended to help assess the capability of thje heterodyne 
radiometer, as presently set up, to yield the quality of spectroscopic data we will later need for the 
design and analysis of upper atmosphere and astronomy experiments. 

To decide this question we will try to extract information on line strengths and widths from 
some of the spectroscopic data so far obtained, then compare this information with known values 
of these quantities for the ob!?erved lines. If the agreement is good enough, we could conclude that 
the heterodyne radiometer is able to measure the kind of spectroscopic data we will need. As will 
become evident from the following discussion, it turns out that the instrument is able to measure 
such data, but improvements are recommended. 

CO -691.473 GHz 

For this line the first study, using AOS I, yielded a set of data which was reasonably well 
calibrated but spanned only the pressure range from about 100 to 300 mtorr. Because the data 
from the third study, using AOS I with an improved data collection and analysis system, had better 
frequency resolution (0.3 MHz as compared to- 1 MHz) and spanned the pressure range from 2 to 
1 600 mtorr, only that data will be used in the analysis given below. The data from the second 
study, using AOS II, was not calibrated and so could not be analyzed. 

The data from the third study yielded a value for the dipole moment of CO which agrees 
quite well with the value given in the literature. The analysis leading to this result follows. 

From the plots of the absorption line at various pressures it is possible to read the line-center 
absorption coefficient and the linewidth at each pressure; the values thus found are given in 
Table 2. 
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Table 2. 

Experimental Results from the Third Study cf the CO Absorption Line at 691.473 GHz 


Pressure 

(mtorr) 

Line-Center Absorption 
Coefficient 
(cm“‘) 

Linewidth 

(MHz;HWHM) 

Total 

Absorption 
(cm"^ MHz) 

2* 

0.00041 

1.9 

0.00078 

5 

0.0010 

2.1 

0.0021 

10 

0.00097 

3.0 

0.0029 

20 

0.0016 

3.2 

0.0051 

50 

0.0022 

3.2 

0.0070 

100 

0.0030 

3,0 

0.0090 

200 

0.0053 

3.6 

0.019 

400 

0.0074 

4.6 

0.034 

800 

0.012 

6.9 

0.083 

1600 

0.015 

12.4 

0.187 


The linewidth is read directly from the plot; the line-center absorption coefficient is cal- 
culated from the ordinate y shown on the plot (labeled “% absorption”) by the relation 



1 

2Q-1 



1 

1-y/lOO 


(19) 


where the first equation is just Eq, 1 (see the Appendix for a derivation) and the second equation 
expresses the absorption coefficient a in terms of the plot ordinate y, which is related to the meas- 
ured quantity Q by y = (1H3) X 200. The length 2 of the sample cell is 190.5 cm. The quantity 
called “total absorption” in Table 2 is the product of the line-center absorption coefficient and the 
linewidth; it is proportional to the area under the absorption line, which is the total absorption. 
This quantity should be proportional to the pressure for all pressures considered here. To better 
see the variation of the line-center absorption coefficient, linewidth and total absorption with 

*This data point wai meaauied in a separate experiment. 
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pressure, these three quantities are plotted vs. pressure in Figure 17. (The common ordinate in | 

Figure 17 is just a relative scale, chosen so that all three quantities fit on the same plot; the three | 

curves were just drawn freehand — no mathematical curve-fitting was done.) i 



Figure 17. Variation of Line-Center Absorption Coefficient, 
Linewidth and Total Absorption with Pressure 
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From Figure 17 It is seen that the linesvidth flattens out at the lower end of the pressure 
range and becomes linear at the high end, while the line-center absorption coefficient rises linearly 
at the low end and flattens out at the higli end and the total absorption varies quite linearly with 
pressure through the entire range. (The two data points at 2 and 5 mtorr are appoi-ently not as 
good as the others, and may bo disregarded.) From the variation of these quantities with pressure 
It is apparent that the lower end of the pressure range approaches the doppler regime and the high 
end is in the pressure-broadened regime. 

From the fact that the data point at 1600 mtorr is in the pressure-broadened regime it follows 
that we may compute the pressure-broadening coefficient Aq from this data point witl; reasonable 
accuracy : we find: 

Aq (300 K) a 7.8 MHz/torr (HWHM) (20) 

where the 300 K figure refers to the fact that the gas in the sample cell is assumed to be at approxi- 
mately room temperature (recall that Aq depends on temperature). With this value of Aq and the 
values of the line-center absorption coefficient given In Table 2 we may use Eq. 9 to compute the 
line strength S at T = 300 K. The only other quantity we need is A^, the doppler halfwidth, which 
is needed to compute the values of x from Eq. 10; A^j is easily computed from Eq. 5, with T = 300 
K and M = 28. We find: 

A^ (300 K) = 0.8 1 MHz (HWHM) (2 1 ) 

If sve regard the doppler region as defined by Aq = pAo < 1/10 Ajj, tlien, using the value of A^ 
found above, we find that the doppler regime holds for p < 10 mtorr. For p = 100 mtorr, A^ a A^ , 
so that for pressures less than about 100 mtoiT the linewidth levels off and decreases slowly until 
Aq « A^ , below which it remains constant. From Table 2 it is evident that the linewidth does 
behave this way (although the numerical values, as might be expected, are not exactly consistent, 
probably due to errors in the measurement of pressure). 
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With tlie values of Aj and Aq computed above we may use Fq. 10 to compute the x-values. 

Eq. 9 may then be used to compute the values of S(300K) from the values of the line-center absorp- 
tion coefficient for the various pressures; the probability integral may bo obtained from tables 
or computed from one of the series expansions given above. The result of these calculations is a 
collection of values for S which range from 77 to 268 Hz cm'"' (dyne/cm^)"A Since S is in- 
dependent of pressure, the approximately 3.5 to 1 variation in these values of S appears to indicate 
experimental errors or shortcomings in the theoretical expressions, or both. In any case the varia- 
tion Is small enovigh to give us an o^de^of-magnitude estimate of S. Wo will therefore simply take 
the mean of the S-vaiues and use it to compute the dipole moment. The mean is found to bo: 

(300K)« 168 Hzcm-‘ (dyne/cin^)"' (22) 

To compute the dipole moment we use Eq. 13, where, for a diatomic molecule such as CO at 
ordinary temperatures, the quantity f, (T) I lUy |- is given by Eq. 17. The resulting relation is: 

S(T) = ( 1 - f-’ (23) 

3c(kT)^ 

where m is the dipole moment of the transition in statcoulomb-cm. If we set! “ 300K, put in the 
values of the mathematical and physical constants (see the Theoretical Background section of this 
report) and express m in debyes (1 debye “ 10“'^ stalcoulomb-cm), we find: 

(1 a 0.07 debye (24) 

Tills value comptires quite well with the value 0. 1 debye given in the literature (see, for example. 
Reference 30). We may thus conclude that the experimental CO data yield a rather accurate value 
for the line strength of the transition. 

HjOa - 601.885 GHz, 606.717 GHz 

For this pair of lines tlic second set of obsei-vations, using .\OS II, yielded line strengths in 
order-of-magnitude agreement with theoretical values (Ref. 36). 

The observations at 5 and 10 millitorr appeared to be in the doppler regime. From the 
observations at 5 millitorr (see Figure 10a), and using the formula 
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( 25 ) 



I 

2Q(0)-I 


(from Eq. I , setting e = 0) to get the line-center absorption coefficient for each line, we found that: 

a(0) = 0.0012 cm"* for the 601 GHz line 

and (26) 

a(0) = 0.0028 cm"* for the 606 GHz line, 

where botli absorption coefficients are for Po = 5 millitorr and Tq = 300 K. If we now use Eq. 3 
with e = 0, p=Po=5 mtorr and T=To=300 K, we can write: 


Sexp (300K) = 


a(0;Po,To) 


Po Sq (O^Jo.Tq) 

for each line. We easily calculate that 

I 0.622 MHz (601 GHz line's 


and 


A.(300K) 


gd(O-.Po.To) = 


SexpOOOK) = 


[ 0.627 MHz (606 GHz line), 

( 0.756 X I0"« Hz"* (601 GHz line) 

I 0.750 X 10"« Hz"' (606 GHz line) 

' 238 cm"* Hz (dyne/cm^)"‘ (601 GHz line) 
,560 cm"* Hz (dyne/cm^ )"* (606 GHz line). 


(27) 


(28) 


(29) 


(30) 


We now want to compare these e.xperimental values of the line strength with the theoretical 
values from Ref. 36; 

[ 0. 1 1 19 cm"^ atm"* (601 GHz line) 

S:^h(300K) =1 (31) 

[0.2045 cm"^ atm"* (606 GHz line) 

To make the comparison we must note that S* and our S do not differ merely in tlie units used, 
but are different quantities. S* is defined in terms of the wavenumber of the radiation, whereas S 
is defined in terms of the frequency. It is not difficult to derive the relation between Che two 
quantities; it turns out to be: 
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S*(T) = 0.338 X 10-^ S (32) 

where the units of S* are cm atm and the units of S are cm“*Hz(dyne/cm-)~^ From Eq, 32 

we may write the theoretical vaiues in terms of S: 

f 33 10 cm Hz (dyne/cm (601 GHz line) 

Sth(300K)= (33) 

(6050 cm“‘ Hz (dyne/cm^)"‘ (606 GHz line) 

Comparison of these values with those given by Eq. 30 shows that the theoretical line strengths and 
the line strengths observed with the heterodyne radiometer are about an order of magnitude apart. 
Considering the necessarily approximate nature of the theoretical calculations for the H^Oj mole- 
cule, which is an asymmetric rotor, and the preliminary nature of the experimental observations, 
these results must be considered reasonably good. 

Sensitivity of Radiometer 

We can determine whether the radiometer is able to detect a small signal approaching the 
theoretical limit by comparing the signal strength of one of the weaker observed spectral lines with 
the minimum observable signal as calculated from the above equation. For this purpose we choose 
the CO line shown in Figure 13a, because at a pressure of 2 mtorr this line is approaching the noise 
level and the good calibration of the plot allows the signal level to be calculated. 

To find AT from Figure 1 3a we must first relate AT to the ordinate y shown in the figure. 

This can be done by using Equations 19, A7 and A9; simple manipulation of these three equations 
yields the relationship we want: 

“ JQQ 

where Tgg is the temperature of the blackbody source, T^^g is the temperature of the gas and 
Tj^gP is the temperature of the reference body (the chopper blade), all in degrees Kelvin. The 
factor j3 may be regarded as including losses in the sample cell and optical coupling losses resulting 
from mismatch between the beam incident on the wliisker antenna of the Schottky diode detector 
and the antenna pattern; preliminary measurements indicate that j3 0.8. 
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From Figure 13a we see that y 5 at tlie center of the absorption line; assuming that 
TgB = 1.300 K, = Tj^ep = 300 K and /3 = 0.8, we find that AT = 10 K. Tlius the radio- 
meter in this case was able to clearly see an apparent temperature difference of ten degrees Kel- 
vin. To compare this with the minimum detectable AT we note that an integration time of 160 
seconds was used for the spectrum shown in Figure 13a; if we assume a noise temperature of 
10,000 K and note that the predetection bandwidth is 0.3 MHz we find from the radiometer 
equation given on page 20 that (AT)j^jj^ = 3.2 K. Tlius the radiometer was able in this case to 
clearly see a signal which is only three times the theoretically minimum detectable signal. 
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SUMMARY AND CONCLUSIONS 


The heterodyne radiometer was used to study selected lines of four gases: hydrogen peroxide 
(HjOj), carbon monoxide (CO), ammonia (NH3) and water (HjO). All lines studied were in the 
frequency range 572-691 GHz, which corresponds to a wavelength range of 434-524 micrometers. 
Two gases were used in the submillimeter laser to obtain appropriate local oscillator frequencies: 
methyl fluoride (CH3F) and formic acid (HCOOH). In all, eight lines were investigated. For three 
of these (a CO line and two H^Oj lines) the data were good enough to do some analysis, giving 
good agreement with theory. From the results obtained in these preliminary investigations it was 
concluded that the heterodyne radiometer is capable of being used for precise spectroscopic investi- 
gations in the laboratory in the submillimeter region of the spectrum, for both absorption and 
emission spectra. 

The prospects are equally good for use of the heterodyne radiometer in astronomical investi- 
gations, particularly the detection of molecular species in interstellar nebulas. The radiometer was 
successfully used to detect the presence of CO in the Orion nebula by observing the J=5-»-6 transi- 
tion at 691 GHz; these observations were made at the NASA Infrared Telescope Facility on Mauna 
Kea, Hawaii, at an altitude of 4,200 meters. Detailed accounts of these experiments are given in 
References 1, 37, 38 and 39. 

These initial investigations using our submillimeter laser heterodyne radiometer system to- 
gether with an acoustooptic spectrum analyzer have generated further development toward a 
second generation unit which is more compact and potentially more efficient for a flight oriented 
package. Design information has been produced to enhance the current AOS I and AOS II devices 
into a rack-mounted version with a me?!"; compact data analysis package. 
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APPENDIX 


DERIVATION OF THE RELATION BETWEEN THE OBSERVED 
DATA AND THE ABSORPTION COEFFICIENT 

The spectra shown in this report are actually plots of a quantity Q which is a simple function 
of the (integrated) voltages at the outputs of the spectrum analyzer channels. What we really want, 
of course, is a plot shewing the value of the absorption coefficient for each channel. We therefore 
need to find the relationship between the measured quantity Q and the absorption coefficient. To 
do this it is necessary to follow in detail what happens to the signal as it passes through the radiom- 
eter. It is convenient to start by finding the relation between Q and the apparent temperatures seen 

at the input to the Schottky diode; these will then be related back to the absorption coefficient. 
The measured quantity Q is defined by: 


(S-R)/R 

(So-R)/R 


(Al) 


where S, Sq and R all represent voltages at the output of any given one of the spectrum analyzer 
channels; S denotes the output voltage when the gas sample cell is filled with an absorbing gas (S 
stands for “signal”), So denotes the output voltage when the gas cell is evacuated and R denotes the 
output voltage when only a reference signal is incident on the whisker antenna of the Schottky 
diode. We now need to express S, Sq and R in terms of the apparent temperatures at the input to 
the Schottky diode; this requires tracing the signal through the radiometer, from the input of the 
diode to the output of the spectrum analyzer. To do this we first need to understand clearly the 
double-sideband mixing action of the diode. 


The output current of the diode is a highly nonlinear function of the input voltage. For any 
nonlinear device, the output corresponding to an input which is the sum of two signals consists of 
products of powers of the two signals, as can be seen by expanding the output as a Taylor series in 
the input variable. If the simple product of the two signals is dominant in the output (except 
possibly for a constant term and the signals themselves), the device can be used as a multiplier. If 
One of the two signals is a sinusoid, the effect of the device is then to merely shift the spectrum of 


Al 
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the other signal both up and down along the frequency axis, by the frequency of the sinusoid. 

Used in this manner, the device functions as a mixer; i.e., as a downward frequency shifter. The 
other signal is then available in a lower frequency range, but with its spectrum otherv/ise unchanged. 
The situation for our device is indicated in Figure Al. 



Figure Al. Input Spectrum to the Mixer 

In Figure A 1 the horizontal axis shows frequency increasing toward the right. The subscripts c, 

LO and IF denote the absorption line center frequency, the local oscillator frequency and the first 
IF frequency, respectively. The wavy line indicates the blackbody spectral distribution, from what- 
ever blackbody source is being used, and the dip in the spectrum at the left represents the absorp- 
tion line being investigated. The local oscillator frequency is shown as being above the absorption 
line, although it may in some cases be below it; By the Raleigh-Jeans approximation to the black- 
body spectral distribution (which is valid for submillimeter wavelengths), 

B = 2kTA" (A2) 

where B is the brightness, k is Boltzmann’s constant, T is the absolute temperature of the blackbody 
source and X is the wavelength. Because the percentage variation in the wavelength is very small 
over the frequency band indicated in Figure A 1 , it is clear that B is nearly constant over the fre- 
quency band. 

The action of the Schottky diode mixer can be imagined as shifting the entire spectrum shown 
in Figure A I down along the frequency axis so that the local oscillator frequency f^o is at zero 
frequency. If the local oscillator frequency lies above the absorption line, the spectrum of the line 
can be imagined to reflect off the zero-frequency vertical axis (if we consider the spectrum to be 
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mathematically represented by only positive frequencies) and to move upward again until the line 
center comes to rest at the frequency ^IF “ ^LO (the first intermediate frequency). The absorp- 
tion line spectrum then overlays a pure noise spectrum (i.e., the blackbody source spectrum) 
in the band fjp ± A. This situation is depicted in Figure A2. If the local oscillator frequency 
happens to lie below the absoiption line the situation is essentially the same, except that in this 
case the line spectrum is not transformed into a mirror image of itself. If we assume that the 
absorption lines we will be dealing with are symmetrical, then we can ignore any mirror inversion 
of the spectrum. 



IF 

Figure A2. Down-Shifted Spectrum 

This intermediate frequency spectrum is then amplified by a low-noise amplifier and passed 
through a filter which passes only the band fiF - ^ (shown darkly shaded in Figure A2). It is 
important to note that this band contains two overlapping spectra: the absorption line spectrum 
plus pure noise. In this sense the action may be described (loosely) as double-sideband detection. 

The band fjj ± A is then down-shifted a second time, but this time it is not altered (except 
for a scale factor), because the noise outside of the band has been filtered out. The resulting band 
± A is then passed through a spfctrum analyzer which perfo.nms a discrete multi-channel spectral 
analysis of the band (and thus of the ab..orption line within it). 

We are now in a position to relate the voltage of the output of the spectrum analyzer to the 
apparent temperature at the input to the Schottky diode. The voltage at the output of any given 
channel of the spectrum analyzer comes from two sources: thermal (blackbody) radiation at the 
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input to the Schottky diode and noise generated within the receiver itself, from the diode on. 
(Note that the thermal radiation is itself just noise; the spectral line “signal” is just a variation 
of the amplitude of this thermal noise along the frequency axis.) If we denote the voltage at the 
output of a spectrum analyzer channel by v(e), where e is the frequency deviation from the center 
of the spectral line, we may write: 

v(e) = h T(e) (A3) 

where T(e) is the apparent temperature at the input to the diode and h is a factor which takes 
account of all gains and losses through the receiver, from the diode on, including conversion losses 
in the diode. Because the two “sidebands” (i.e., the band containing the absorption line and its 
image on the other side of the local oscillator frequency) differ by only a small frequency, we may 
assume that they are affected to the same degree by conversion losses in the Schottky diode; at 
all points downstream in the receiver the two sidebands are in the same frequency band, so that 
they undergo the same losses and gains. We may thus assume that h does not depend on frequency 
over the band we are working with; note also that h includes dimensioned conversion factors, since 
V and T have different dimensions. The temperature T is the sum of two fictitiou’; wmperatures: 

• The temperature of a fictitious 'blackbody which, if viewed through empty space, 
would give the brightness actually incident on the Schottky diode; this brightness 
is related to T by the Raleigh-Jeans Law (Eq. A2). 

• The temperature of a fictitious blackbody which, if its radiation were incident 
(through empty space) on the diode, would give rise to the output voltage observed 
in the absence of any true input signal (i.e., the voltage caused by internally generated 
noise). In other words, this component of T is just the receiver noise referred to the in- 
put and expressed as an effective blackbody temperature. 

We now must write explicit forms for T, for three distinct cases: 
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i) The gas sample cell is filled, so that an absorption line is present; 

ii) The cell is evacuated ; 

iii) The input to the diode comes from a reference source. 

For tlie first case we may represent the part of T due to the brightness incident on the diode by 
*^OFF"LiNE ^"^OFF*LiNE where Tqpp.[^ij>i£ denotes the apparent temperature seen by 

the diode away from the absorption line and AT(e) denotes the difference in apparent temperature 
between the off-line level and a point on the line; note that Topp.^ipiE appears twice because of 
the two overlapping sidebands. As will be described in more detail below, Toff-line the actual 
temperature of the blackbody source, reduced by internal losses in the gas ceil walls and associated 
windows and mirrors; AT(e) is the temperature difference due to absorption, also reduced by the 
same internal loss factor. The part of T due to internally generated noise (i.e., the noise tempera- 
ture of the receiver) may be represented by 2T(^, since tlie noise temperature may be assumed to 
be constant over the frequency band considered; again the factor of two is due to the two side- 
bands. Thus we may write T for the first case in the form: 

T=‘2Topp.uNE + 2TN-AT(e) (A4) 

The second case is the same, except that AT is zero: 

T = 2Toff,line 2T^ (A5) 

For the third case we have: 

T=2Tapp+2TN (A6) 

where Tref is the actual physical temperature of a reference body placed riglit before the receiver 
(actually, before the diplexer). If we put the tliree expressions forT given by Eq’s. A4, A5 and A6 
into Eq. A3, denote the resulting v’s by S, So and R, respectively and put those into Eq. Al, we get 
the measured quantity Q in terms of the apparent temperatures seen by the receiver. Before writing 
out this relation, however, it is necessary to briefly describe the measurement process as it is carried 
out in the laboratory. A chopper (a fanlike device with two or more blades) is inserted into the 
beam just before the diplexer. The rotating blades interrupt tlie beam emerging from the sample 
cell several times a second, so that the receiver alternately sees the beam from the cell and one of 
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the chopper blades, which are at room temperature; these blades supply the reference, Tref* 
Observations are first made with the sample cell filled with the absorbing gas, for a long enough 
time to average out the random fluctuations and obtain a good spectrum. The sample cell is then 
evacuated and the process is repeated. If the elapsed time between the two sets of observations is 
too great, the value of h (in Eq. A3) may have changed due to drifts in the gains and losses through 
the receiver. <We assumed that h is independent of frequency over the small range we are working 
with, but it may vary with time.) Because of this possibility we will attach subscripts to n and to 
Tpi ; the subscript 1 refers to the first set r f observations with gas in the cell and the subscript 2 
refers to the second set of observations with the cell evacuated. Thus we may write: 

S = hi(2ToFF.LiNE ” AT(e)) 

Sq ^2(^"^0FF.LINE ) 

^ 1^) 

where the subscripts in the expression for R depend on whether R is measured in the first or second 
set of observations. If we put these expressions into Eq. A1 we obtain the relation between the 
measured quantity Q and the apparent temperatures seen by the radiometer; before simplification 
this relation takes the form; 

th,(2TopP.UNE +2T„. - AT(e))-h,(2TReF . 2T>,,)l/h,(2T^Er + 

[hj (2TofF-LINE ■*'2Tj^2) “^2(2^REF /^2(2^REF "^^T^j) 


From this formula for Q it is apparent why Q was defined the way it was: it is so that all the h’s and 
Tj^’s cancel, provided (and this is essential) that the elapsed time between the observations with gas 
in the cell and with the cell evacuated is small enough so that no appreciable drift has occurred in 
the receiver between the two sets of observations. If this is the case then the quantity Q is inde- 
pendent of the noise temperature of the receiver (and of the gain or loss in any component). The 
expression for Q then takes the simplified form: 


Q = 


1-/2 


AT(e) 


("^OFF-LINE "^REf) 


(A7) 
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I t Is important to note one other property of the definition of Q: the value of Q remains the same 
under any linear transformation of the output voltage v. Thus if Eq. A3 were to bo replaced by a 
relation of the form v(e) » aT(e) + b, the same final formula for Q would have resulted. This 
means that the value of Q does not depend on the zero ref.-rence used to define the ouput voltage 
V, or on any change of scale factor, provided only that v is defined the same way in measuring all 
three of tire quantities S, So and R. 

We now must relate the quantities appearing in Eq. A7 to tlie absorption coefficient of the gas 
being investigated. To do this we need to apply the theory of radiation transfer to the gas cell with 
tlie blackbody source at one end, to calculate tlie apparent temperature at the other end Oust 
before the diple.\er). The theory (too lengthy to give here) yields the following relation: 

T ^ «LOSS)fi + T.. . <^(1 ~ + «LOS5>S (ASl 

where a(e) is th*? absorption coefficient, a^oss coefficient due to internal losses in the 

gas cell (wall absorption plus losses in windows, etc.),fi is the length of the cell, T(e) is the apparent 
temperature at the far end of the cell (just before the receiver), Tgj, is the physical blackbody 
temperature and Tq,^s is the physical temperature of the gas in the cell. The first tenn on the 
riglit-hand side of Eq. A8 is due to absorption by the gas; the second term takes account of emission 
from the gas, assuming local thermodynamic equilibrium. The temperature of the blackbody source 
was purposely made quite high, so that the absorption tern\ would be dominant and the emission 
term could be neglected. However, we will carry the emission term througli the calculations, for 
the sake of clarity. Using Eq. AS we may immediately write e.xpressions for the off-line apparent 
temperature and the apparent line depth: 

Toff.une =T„„ e-“tossS +To.,s(1 -e-">-oss») (A9) 

and 
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AT(e) =*Toff.line “'Trec(^) 

(A 10) 

= e““Loss£ (Tbb - To AsX I - ) 

If we put the two expressions given by Eq’s. A9 and AlO into Eq. A7, and if we write (1 “j3) = 
e~*Loss * (so that j3 is just the fractional loss of apparent temperature in the gas cell due to internal 
losses), we obtain: 

Q(e)=l-V^ = (All) 

1 ^GAS ^REF 

(1-13)(Tbb-Tgas) 

It is clear that the fraction in the denominator must be small, since Tp, £p is about 300 K (room 
temperature), Tgg is about 1274 K and Tq^s is somewhat above room temperature due to heating 
of the gas (/3 = 0,4); we may thus write, to a good approximation: 

Q(e)2s'A(l +e-“<«>2) (A 1 2) 

Note that we could have arrived at this same approximation by neglecting the second term (the 
emission term) in Eq. AS. Also it is worth nothing that Vi<Q^ 1, because 0 < a. If we solve 
Eq. A 12 for tlie absorption coefficient we obtain the relation we wanted: 


a(e) =*~ln 

. 2 2Q(e)”l 


(A 13) 


This gives the absorption coefficient in terms of the measured quantity Q. If the values of the 
absorption coefficient thus obtained are plotted for each of the spectrum analyzer channels, the 
result is the absorption line shape, to within the resolution of the analyzer. 



